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Abstract 

We calculate the second virial coefficient of spin-1/2 anyon gas in the various 
values of the self-adjoint extension parameter by incorporating the self-adjoint 
extension method into the Green's function formalism. Especially, the com- 
pletely different cusp- and discontinuity-structures from the result of previous 
literature are obtained when the self-adjoint extension parameter goes to in- 
finity. This is originated from the different condition for the occurrence of 
irregular states. 
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Since the Kronig-Penny model |T[ described successfully the band structure of energy 
spectrum in the solid-state physics, the point interaction problem has been applied to the 
various branches of physics for a long time. Recent application of the point interaction in 
the theoretical physics seems to be concentrated upon the subjects related to the quantum 
mechanical renormalization and anyonic theories. [Q 

The frequent use of the point interaction in the quantum mechanical renormalization is 
mainly due to its advantage of permitting the derivation of the exact solution. Hence, the 
comparision of an exact solution with a perturbative solution, which can be obtained by solv- 
ing the Lippmann-Schwinger equation iteratively, allows us to understand the subtleties of 
renormalization scheme encountered frequently in the quantum field theories. It may be also 
very helpful to understand the highly non-trivial concepts like dimensional transmutation 
and anomalies 0. 

In spin- 1/2 Aharonov-Bohm problem which is directly related to the fermion-based 
anyonic theory, the two-dimensional point interaction is realized as a Zeeman interaction 
of the spin with a magnetic flux tube. The two different approaches, renormalization and 
self-adjoint extension [0], for the quantum mechanical point interaction and the coincidence 



of their results are presented by Jackiw |TU[. In Ref. W^ one of us analyzed the incorporation 
of the self-adjoint extension into the Green's function formalism by solving the Lippmann- 
Schwinger equation nonperturbatively. 

One of the way to investigate the statistical properties of the spin-1/2 anyon system is 
to evaluate the second virial coefficient as a function of a statistical parameter. The second 



virial coefficient of the Boson-based anyonic system is firstly calculated by Arovas et al. [[12 
The most interesting feature of their result is that the virial coefficient interpolates between 
the Bose and Fermi values with a periodic dependence on the flux a carried by anyon. 
Subsequently, the second virial coefficient for a system of spin-1/2 anyon gas is computed 
by Blum et al. [|13| by using the condition of the irregular solution 

I m + a |< 1, 
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I m I + I m + a 1= —as, (1) 

where m is the angular momentum quantum number and s is ±1 for spin up and spin down, 
respectively. Their striking result is that there exist discontinuities in the virial coefficient 
at all even, nonzero values of a. 

In this Letter we compute the second virial coefficient for a spin-1/2 anyon gas by in- 
corporating the self-adjoint extension method into the Green's function formalism. In this 
case, it is well-known that the arbitrary combination of the regular and irregular solution, 
which is characterized by the self-adjoint extension parameter, is derived when | m + a |< 1. 
We expect the absence of the latter condition in Eq.(|l|) for the appearence of the irregu- 
lar solution in the self-adjoint extension method does yield completely different cusp- and 
discontinuity-structure in the virial coefficient from the known in the previous literatures. 
The evaluation is performed at various self-adjoint extension parameter X^- In particular, 
we find new discontinuities at all non-zero integer values of flux when —>■ oo. When, 
however. Am becomes finite, the discontinuities turn into cusps. This might be an existence 
of different kind of phase transitions from the previous cases. 

In general, the second virial coefficient can be calculated by making an energy spectrum 
discrete. Therefore, we introduce the simple harmonic oscillator potential in the Hamitonian 

H = Ho + v6{r), (2) 

where 

The energy-dependent Green's function G'[ri,r2 : E] for H can be derived by solving the 
Lippmann-Schwinger equation nonperturbatively. Following Ref. [|ll|] one can arrive at 



iKr.:B] = ^j!^l^^^te^. (4) 

- + lim,2^,^+ Go[e2,ei : E\ 

where 

i[r2,ri : E] = G'[r2,ri : E] - Go[r2,ri : E] (5) 



and Go[r2,ri : E] is energy-dependent Green's function for Hq [jl^. Since (j[r2,ri : E] 



J2n 4'ni^2)4'n{^i) / {E + En) , G and (p niust obey the same boundary condition, which is given 
by self-adjoint extension method at the origin |T^. In Ref. |jll]] bound state energy is derived 



by imposing the boundary condition to A[r2,ri : E] instead of (j[r2,ri : E]. In this case 
the bound state spectrum can be obtained by solving the equation 

^ ^ r(i+|m + a|)r( ^"i'"%"i+^/" )' ^ ^ 

when I m+a \< 1. Here, is the self-adjoint extension parameter. If Go[^2, ri : E] does not 
have any pole, the imposition of the boundary condition to A[r2, ri : E] is physically relevant 
because the most contribution to (j[r2,ri : E] is A[r2,ri : E] at bound state energies. If, 
however, Golr^,^^ : E] does have poles, we can not conclude a priori the relevance of the 
above procedure. In order to get the credibility in this case, the poles in Go[r2, ri : E] must 
not contributed to the poles in G'[r2,ri : E] at the final stage. We confirmed by explicit 
calculation that this is indeed the case in our model. Furthermore, imposing the boundary 
condition to G'[r2,ri : E], one arrives at the same result. 

When I m + a |> 1, no irregular solution occurs, and therefore, the bound state energies 
are given by 

en,m = (,'^n + 1+ \ m + a (7) 

where n is a non-negative integer and m is an integer. We first consider the two special cases 
in Eq.d): A^ = and A^ = oo. In both cases, we can obtain the analytic solutions for the 
bound state spectrum. For the former case, the bound state energies coincide with Eq.(0), 
which means that there is no irregular states even at | m + a \< 1. (hard-core repulsion 
potential exists.) For the latter case, however, the bound state energies are given by 

e„_m = (2n + 1- I m + a (8) 



This is exactly the same as the case of Blum et al. The only difference is that there 



is no I m I + I m + a 1= —as condition for the occurrence of the irregular solution in our 



self-adjoint extension method. We plot the energy eo,o for Am = cc as a function of a in Fig. 
|I|. In contrast to Ref. |]T3|, there are two discontinuities at a = —1 and +1, respectively, 
and one cusp at a = 0. This will probably imply the different behavior of second virial 
coefficient from that of Blum et al. 

When < Xm < oo, the bound sate spectrum for | m + a |< 1 can be found by graphical 
analysis. The splitting between the adj ascent levels is not 2uj unlike the two above-mentioned 
cases, but depends on Am and a. Using the energy eigenvalues obtained numerically, we now 
calculate the second virial coefficient 



.(a=0) 



(9) 



where B^iT) is the second virial coefficient of free Fermion (Boson) which is given by \}? 
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Xj^), At is the thermal de Broglie wavelength and /? = 1/kT. 



First let us consider the Boson-based case. We perform the summation only over the 
even m, and obtain 



B^{T)-B','^{T) 
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even, 



= odd, 



(10) 



^ En,|m+a|<l 

where a = + 5, is an integer and < 5 < 1. In the a; — limit, the result for 
i?^(T)/A|. at different self-adjoint extension parameters are shown in Fig. ^ The simple 
numerical calculations are needed to perform the summations when | m + a |< 1. We add 
the e"'^""-'" factors up to n = 5, for which the errors are less than 1 %. When Am = 0, the 
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result is exactly same as that of Arovas et al. When Am begins to be finite, the cusps 



at the Fermion points begin to appear and become more deeper with increasing Am- This 
fact agrees with the recent calculation of second virial coefficient of anyons without hardcore 
p!6| . The appearence of cusps at the Fermion point is interpreted in Ref. as follows: 
as the possibility of the overlap between particles is introduced, the particles becomes more 
Boson-like. The cusps at the Bose points still remain, at variance with the result of Ref. 
| 16| which is unbelievable on account of the abrupt appearence of the inflection point at 



e = 10. As the self-adjoint extension parameters increase, the values at all points except 
the Bose point become smaller than the hard-core values. This can be considered as follows. 
The introduction of the irregular states at the origin increases the possibility of overlap 
between particles, so that it makes the system compressed. Increasing the portion of the 
irregular solution in wave function, the augmentation of the self-adjoint extension parameter 
decreases the second virial coefficient. 

The Fermion-based case can be calculated by summing odd m's. We consider only 
the unpolarized Fermion system, which is obtained by averaging over the four spin states 
(triplet + singlet). The virial coefficient -Bf(T)/A|. in the ^ is plotted in Fig. ^ for 
\m = 0, 1/5, 1 and oo. When A^ = 0, the result without irregular solution in Ref. pTS 



IS 



recovered. But for finite Am, the value at Bose point jumps down suddenly to the negative 
value and the second virial coefficient decreases with increasing Am as in the Boson-based 



case. As Am — ^ c>o (the case of Blum et al [T^), there exist two kinds of discontinuity at 



both Boson and Fermion points, which shows the different behavior from Ref. |T^. This 
means that it is some different phase transitions from those of Blum et al. The origin of 
this difference comes from the different condition for the occurence of irregular states. As 
Am decreases, discontinuities turn into cusps at integer points and as Am — > 0, the result is 
reduced to the case in which the irregular solution is ignored. 

In conclusion, using the self-adjoint extension method and imposing the boundary con- 
ditions at the origin, we calculate the second virial coefficients for the Boson-based and 
Fermion-based anyon systems. Both cases show the decreasing behaviors when the self- 
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adjoint extension parameter Am increases. This tells us that the values of parameters de- 
termine the magnitude of repulsion between two particles. Although, for the Boson-based 
case, the decreasing behavior also shows up in Ref. |T6[, the cusp structure at Bose point 
is completely different. It is also shown that the Fermion-based case has a different kind 
of discontinuity from that of Blum et al. |TB[, which is due to the absence of the second 
condtition in Eq. (|l]). This might exhibit another kind of phase transition. 
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FIGURES 

FIG. 1. The bound state energies eo,o as a function of a for = oo case. 

FIG. 2. B^(r)/A|n for Boson-based anyons at various self-adjoint extension parameters A^. 

FIG. 3. B2{T)/\j, for unpolarized Fermion-based anyons at various self-adjoint extension 
parameters A^^. 
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